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Abstract 

Derrida et al. and Schiitz and Stinchcombe gave algebraic formu- 
las for the correlation functions of the partially asymmetric simple 
exclusion process. Here we give a fairly general recipe of how to get 
these formulas and extend them to the whole time evolution (starting 
from the generator of the process), for a certain class of interacting 
systems. We then analyze the algebraic relations obtained to show 
that the matrix approach does not work with some models such as 
the voter and the contact processes. 
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1 Introduction 

A few years ago Derrida et al.EI'El suggested an intriguing "matrix ap- 
proach" to the one-dimensional Asymmetric Simple Exclusion Process (ASEP). 
This approach has later been used to treat variants of the modelEl'I3>IHl^ 
extended to non steady states and by Schiitz et al.^l'^l'^l, used to 
study fluctuations by Derrida et al.EI and the multispecies case by (among 
others) Isaev et al. . 

The main aim of this paper is to study the difficulties that arise in 
potential applications of the matrix approach to cases in which the near- 
est neighbor interaction or the particle conservation (both present in the 
ASEP) are violated. Further light on the applicability of the matrix method 
is shed by the integrability criterion illustrated by Popkov et al.^l. 

In section |5] we provide a general recipe (using the generator of the 
process) to find the algebra of the matrix formalism associated to both the 
steady state and the whole dynamics of any one-dimensional interacting 
system such that at each step the configuration changes only in two adjacent 
sites. A more complete description, with a pedagogical aim will be given 
elsewherel^ . In section|31wc apply the recipe to some important interacting 
systems such as the contact and voter models and show that the matrix 
algebra obtained is not useful to treat them. 

We will consider only systems in the lattice {!,..., N}, this is an in- 
trinsic limitation of the matrix approach. The dynamics of an interacting 
particle system is usually defined by giving the generator of the process, 
the general form of which can be found for instance in Liggett 's bookCSI. 

For example, the generator of the ASEP, if particles jump one site to 
the right (left) with rate p {q = 1 — p) and enter the lattice from the left 
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(right) at rate a (5) and leave it at rate 7 (/3), is defined by 
W)(t) = 

J2 [pr{x){l-r{x + l))+qT{x + 1)(1- r(a:))][/(T-'-+i) -/(r)] 

x=l 

+ [a(l-r(l))+7r(l)][/(ri)-/(r)] 

+ [<5(l-r(iV))+/3r(7V)][/(r^)-/(r)] (1) 

where r = {'''(a^)}^^i is the configuration of the system, r^'^ is the con- 
figuration obtained fi'om r by exchanging the content of the sites x and y, 
and is the configuration obtained from r by changing the content of the 
x-th site. 

In the following formulas, \V)) is a vector in an (as yet) unspecified 
linear space equipped with an inner product, D and E linear operators on 
the same space, {{W\ is an element of the dual space. So is the 

inner product generally written as (W,ylV). 

The formula of Derrida et al. to write the probability of a given config- 
uration in the stationary state of the ASEP isHi 

1 ^ 

where D, E, \V)), {{W\ are matrices and vectors that satisfy 

{(3D-6E)\V)) - \V)) , 

pDE - qED = D + E , (3) 
{{W\{aE~jD) ^ {{W\ 

and Zjv is a normalization factor. 

One can check these formulas provide a sufficient condition for the mea- 
sure to be stationary by observing they satisfy the recursion relations for 
the probabilities (first due to Liggett^^) that relate the probabilities for 
the system with K sites to the ones for the system with K — 1 sitesD. 
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2 From the Generator to Matrix 
Product States 



Let us start by re-writing the generator by making use of a formalism bor- 
rowed from quantum mechanics. For all j = 1, iV let us define the Hilbert 
space Hj := span^Q)- , — . Consider the operators ,n,m 

defined by: a+|0) = |f),a-|0) = ,n|0) = , m I - n ,a+|l) 
,a^|l) ~ |0) ,n\l) = |f), where I is the identity. Interpreting |0} and 
|1) as empty site and occupied site respectively, the role of a+, a^, n as cre- 
ation, annihilation, number operators respectively is rather obvious. The 
most immediate choice of an explicit expression for the operators and vec- 



tors above is |0) 
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0\ 0\ ^jv 

, m = . Now we take the tensor product H^r = (X) _i Ti-j 

1 oj \q ij '~ 

to describe the system on all the N sites. 

If we consider for example the ASEP, in this "quantum hamiltonian" 



formalism 
H 



the generator is given by 

Xlp("fc «fc+i - nkmk+i) + g(a^afc^i - mknk+i)+ 

k 

7(a;f - ni) + a{a'^ - mi) + /3(a^ - njv) + 5{a'^ - ruN) 
^hf + Y,hk + h%, (4) 

k 

where the superscript d denotes a boundary term and 

/o o\ 



7 

-7 



hk = 




p -q 
-p g 
\0 OJ 

For any given operator or vector b in the space Tik we use the notation 
bk = I(8)---(8)I®6(8)I---(8)I with b as k-th factor. Using a different h, , 
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this formulation can be used for any process (Uke the voter and contact, 
e.g.) such that the occupation number of each site is either or 1, and 
such that the dynamics involves a couple of neighboring sites at a time 
(slight generalizations can be treated as wellEIIHI). The generator in is 
the same as in as can be checked by computing the Dirichlet Form for 
both and verifying that they coincide (the same holds for processes with 
different h). It is however easier to look closely at each part and see what 
it does. For instance a'^a~ represents a jump to the right and nm takes 
into account the complementary event (the particle stays where it is). 
We now look for a stationary solution of the master equation 

\Pit)) = H\Pit)) (5) 

which describes the dynamics of the system by giving the time evolution of 
the vector of probabilities of configurations, i.e. we look for a distribution 
\P,) such that H\Ps) = 0. 

In order to show where the general idea can be guessed from, let us con- 
sider again the case of the ASEP, to show that under special conditions 
(namely (a + /3 + 7 + S){a(3 — -f5)/{a + d){f3 + j) — p — q) the stationary 
state is a product state: jPg) = (where d = (a + 5)/{af3 — ^S) 

, e — {P + 7)/(ck/3 — 7(5) and the normalization constant is clearly Zj^ ~ 
{e + d)^). To prove that H\Ps) = 0, one should first check 

This makes the sum through which H is defined telescopic (recall that we 
are omitting the factors of the tensor product on which the operators act 
trivially as the identity), and since 

the cancellation of the first term with the last is assured by the boundary 
terms. 
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In other words, H\Ps) = would be solved for instance if we had zero 
for all i in the r.h.s. of but this is too restrictive, so we look for the 
first non trivial possibility: instead of zero, we impose a "telescopic term" . 
This is inspired by the dynamics, that acts with the same /i. on all couple 
of adjacent sites, so the generator acts twice on each site. We will now 
try to make the above approach work for non-product states by imposing 
a similar telescopic property. The idea is to move into a richer context, 
substituting the numbers 1, e, d appearing in 10 with some time-dependent 
operators (non commuting and acting on an auxiliary space of generally 
infinite dimension) S, E, D to be determined, aiming to get the weights of 
each possible configuration through a bracket with a couple of vectors {{W\ 
and \V)) to be introduced in the same space. For instance for a system 
consisting of a single site we would impose ((W^|(^)|^)) — {\(w\e\v])) ~ 
(^) and clearly, in the case of a product measure, {{W\ {^)^^ \ V)) — (^)'^ . 
We can also write H{{W\Q^^\V)) = 

Let us now write \P) — \V)) for the probability vector 

and plug it into the master equation Clearly Zn — {{W\C'^\V)), with 
C = D + E, that does not depends on time by conservation of probability. 

It is easy to show that the master equation Q is satisfied if the following 
equalities hold (thanks to the same telescopic cancellation mechanism we 
used for the product state) 

Id \ fD\ fD 

+h 



(7) 



2dt. jyEj yE 
e) ^ ts) - i's) « (e 



{{W\ 



Id , fD\ f-S 



^2jt^^^'\E) \S 



These are the relations of the matrix algebra of the process. If we chose for 
example the h, of the ASEP, these equations take the explicit form of the 
algebra found by Stinchcombe and Schiitz^l'^l that includes as a special 
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case the stationary one ^ of Derrida et al. (taking S' = I and putting all 
the time derivatives equal to zero). With this procedure we can exhibit 
an algebra for all the models with a dynamics involving only a couple of 
neighboring sites at a (see the same works for a classification 

of the models with different ft..). If one found an explicit expression for all 
the operators and vectors, the model could in principle be solved exactly 
(provided the algebra is not empty). Unfortunately, this is in general very 
difficult to accomplish (a purely algebraic treatment can also be used^l). 
In the case of the ASEP, thanks to the preservation of the number of 
particles in the bulk dynamics, the local generator h, has a block form, 
with zero entries in the first and last row and column. This special form 
of /i. is such that in stationary conditions the four equations collapse 
to just one: Q- But this great simplification may not occur for different 
models. In many cases the algebra can be empty (or too complicated to 
deal with), as we are going to show for the contact and voter models. We 
can say that the method works for the processes, such as the ASEP, the 
probability measures of which are either product, or a generalization that 
we can classify as "matrix product measures" . If one distinguished only 
between product and non-product states, the choice would be in general 
only between a numerical tensor product and a convex combination of as 
many such products as the cardinality of the configuration space. If the 
states of a process are matrix product, one can chose to deal again with a 
single tensor product, thanks to the richer nature of the entries, matrices 
instead of numbers. 

Algebras defined by conditions like lO are called Diffusion Algebras 
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3 The matrix Approach Beyond 
Simple Exclusion 



3.1 Exclusion process with double jumps 



The method to write the matrix algebra of the process can also be extended 
to the case of dynamics not limited to neighboring sites, such as for instance 
the exclusion process with jumps of length two permitted. The generator, 
in the case of symmetric dynamics, is (up to boundary terms): 



H = 



k 



i"'k ~ nkuik+i) + («fe a^+i - TOfenfc+i)+ 
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For this system we impose the telescopic property to solve the master 
equation in the following way: 



1 d 



-S' 

s 



D 

E 
D 
E 



-S 

s 

-s 
s 



which is the same as 
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-DSD + 2SD^ - D^S 



-{2DDE + DDE + D'^E) + D'^E - ED'^ = -DSE + 2SDE + D'^S 
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-{2D ED + DED + DED) + DED - ED^ = 
3 

DSD + 2SED - DES 

-{2DE^ + DEE + DEE) + 2DE^ - EDE - E^D = 
3 

DSE + 2SE^ + DES 

\{2ED^ + EDD + EDD) - D^E - DED + 2ED^ = 

-ESD - 2SD^ - EDS 

-{2EDE + EDE + EDE) - DE'^ - EDE = 
3 

-ESE - 2SDE + EDS 
\{2EED + EED + E^D) - DE^ + E^D = ESD - 2SED - E^S 

-{2EE^ + EEE + E^E) - = ESE - 2SE^ + E^S 
3 

These relations define now a cubic algebra, as opposed to a quadratic one, 
which is therefore not a Diffusion Algebra in the sense of ^ . Unfortunately 
algebras of degree higher than two are very difficult to handle (see e.g. 
VershikEl). However algebras of degree higher than two appear e.g. in 



3.2 Voter and Contact Models 

For a description of the voter and contact models see Liggett^l . It is easy 
to see that the local generator for the voter model can be written in the 
form of the r.h.s. of I^J with 



/o 






-1 -1 
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-1 



9 




where A and fi arc the rates for opinion changing in the boundary sites. 
Notice that there are non zero entries in the first and last row. It is easy 
to compute that 

/ -{D,E} \ 




2DE 

2ED 



V -{D,E} ) 

and so we can conclude that the algebra and its stationary limit are given 
by 

^{DD + Db) - {D, E} = [S, D] {D, E] = 
]^{bE + DE) + 2DE = SE + DS — > 2DE = C 
^{ED + ED) + 2ED = -{SD + ES) — > 2ED = -C 

1{EE + EE) - {D,E} = [E,S] — > {D,E} = 

Hence in stationary conditions 

[D,E]=C = D + E, {D, E} = 0, iiD\V) = \V), (W\XE = {W\. 

Notice that the relations are similar to the ones of the ASEP, but there is 
an additional condition: D and E anticommute. 
The local generator of the contact model is 



h. = 





—a 


—a 


o\ 
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h 



so that 

/ -a{D,E} \ 
(a + 0)DE - aE^ 
{a + (3)ED - aE^ 
\ -0{D,E} + 2aE^ J 

and so we can conclude that the algebra is given by 



-{DD + DD)-a{D,E} = [S,D] 
^{DE + DE) + {a + (3)DE - aE^ = SE + DS 
^{ED + ED) + {a + P)ED - aE^ = -{SD + ES) 
^{EE + EE) - I3{D, E} + 2aE'^ = [E, S] 

so that in stationary conditions E^ = 0, [D,E] = C ,{D,E] = if we 
assume a = /? = 1. 

Clearly these relations define a subalgebra of the one for the voter 
model. 



Theorem 3.1 In stationary conditions, the algebra of the voter model is 
empty (and a fortiori so is the one of the contact process and so are the 
ones for the whole time evolution). 

Proof The algebra is 

DE={D + E)/2 , ED = -{D + E)/2 , 
DE = -ED , D\V) = fi\V) , {W\E = {W\X . 

If 

'd.=D, E 
we get, from the first two conditions 

N 

{W\ n MV) = {W\[P{D) + Q{E)]\V) = [P(1/m) + Qil/X)]{W\V) 
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with some polynomials P and Q; but the third condition (anticommuta- 
tion) also implies 

N 

{w\ n MV) = {±){w\E"'D^\v) = {±){l/\y^{l/^Jir{w\v) 
fc=i 

where m + n — N . The two expressions cannot be equal for all values of A 
and ^. □ 

This shows that, following the recipe of section ((2Jl, we cannot use the 
matrix approach. However, the l.h.s of |(7J) reflects directly the dynamics of 
the process and does not depend on the matrix formalism, but the telescopic 
r.h.s. is only inspired by the nearest neighbor nature of the dynamics 
and it is more "artificial" . In other words, if another way to solve the 
master equation were developed, some kind of matrix approach could still 
be productive also for those models that cannot be treated with the current 
matrix approach illustrated in this paper. 
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